Pseudospin Anisotropy Classification of Quantum Hall Ferromagnets 
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Broken symmetry ground states with uniform electron density are common in quantum Hall systems 
when two Landau levels simultaneously approach the chemical potential at integer filling factor v. 
The close analogy between these two-dimensional electron system states and conventional itinerant 
electron ferromagnets can be emphasized by using a pseudospin label to distinguish the two Landau 
levels. As in conventional ferromagnets, the evolution of the system's state as external field param- 
eters are varied is expected to be sensitive to the dependence of ground state energy on pseudospin 
orientation. We discuss the predictions of Hartree-Fock theory for the dependence of the sign and 
magnitude of the pseudospin anisotropy energy on the nature of the crossing Landau levels. We 
build up a classification scheme for quantum Hall ferromagnets that applies for single layer and 
bilayer systems with two aligned Landau levels distinguished by any combination of real-spin, orbit- 
radius, or growth direction degree-of-freedom quantum numbers. The possibility of in-situ tuning 
between easy-axis and easy-plane quantum Hall ferromagnets is discussed for biased bilayer systems 
with total filling factors v — 3 or v = 4. Detailed predictions are made for the bias dependence of 
pseudospin reversal properties in v — 3 bilayer systems. 
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I. INTRODUCTION 

Studies of magnetic phenomena in semiconductors 
have opened fruitful new ways to explore the subtleties 
of quantum magnetism. In quantum Hall samples, the 
tunability of semiconductor electronic systems and the 
quantization of single-particle energies into macroscop- 
ically degenerate Landau levels (LLs) combine to open 
up a rich and varied phenomenology. The effective zero 
width of electronic energy bands enhances the role of 
inter-particle interactions!!] and can frequently lead to the 
formation of ordered many-particle ground states, includ- 
ing ferromagnetic ones. n 

Most studies of quantum Hall ferromagnetsB (QHFs ) 
have focused on spontaneous spin- alignment in a single- 
layer two-dimensional (2D) electron system at LL filling 
factor v = 1. In this case it turns out that electron- 
electron interactions favor fully aligned electron spj. 
even in the limit of vanishingly small Zeeman coupling 
and the ferromagnetic ground state of the system is de- 
scribed exactly by Hartree-Fock (HF) theory. Because of 
the near spin-independence of the Coulomb interaction, 
the v = 1 single layer QHF is a Heisenberg-like isotropic 
two-dimensional ferromagnet. One of the unique prop- 
erties of this simple itinerant electron ferromagnet is 
that its iastantonsa (Skyrmions) carry charge and can be 
observedQ in the ground states at filling factors slightly 
deviating from 1. 

The notion of the QHF can be generalized however. It 
turns out that, at least according to HF theory, broken 
symmetry ground states occur at integer filling factors in 
quantum Hall systems any time two or more valence LLs 



are degenerate and the number of electrons is sufficient 
to fill only some of the LLs. In effect, electrons in the 
ordered state occupy spontaneously generated LLs that 
are linear combinations of the single-particle levels chosen 
to minimize the electron interaction energy. 

The simplest example of a pseudospin QHF obtains at 
v = 1 in balanced bilayer 2D systems where the single- 
particle LLs in the|-t.sKO-layers are degenerate. In the or- 
dered ground state&Era the electrons occupy a LL which 
is a linear combination of the isolated layer levels, form- 
ing a state with spontaneous inter-layec phase coherence. 
Recently, Josephson-like behavior seenEI in 2D-to-2D tun- 
neling spectroscopy studies of bilayer systems has pro- 
vided a direct manifestation of collective behavior gen- 
erated by this broken symmetry. In v = 1 bilayer sys- 
tems, the broken symmetry state minimizes the Hartree 
energy cost by distributing charge equally between lay- 
ers and gives up part of the intra-layer exchange energy 
while gaining more in inter-layer exchange energy. Un- 
like the single layer v — 1 QHF, the ordered HF ground 
state is not exact in this case, and the order predicted 
by HF theory can be destroyed. With decreasing inter- 
layer exchange energy quantum fluctuations around the 
mean-field ordered state become more important and for 
layer separations larger than approximately two magnetic 
lengths fluctuations destroy the spontaneous coherence. 
The corresponding order-disorder quantum phase transi- 
tion has been observed experimentallyu In v = 1 bilayer 
QHFs, it is the layer degree of freedom which is repre- 
sented as a pseudospin-1/2. With this mapping the phase 
coherent state is equivalent to a spin- 1/2 easy-plane fer- 
romagnet. Finite temperature Kosterlitz-Thouless phase 
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transition!, continuous quantum phase transition in- 
duced by in-plane magnetic fielda, and macroscopic col- 
lective transport effectso are are among the remarkable 
phenomena which have been studied|-on bilayer QHEs. 

Recent experiments in single layerl!3 and bilayertil 2D 
systems at even-integer filling factors have further en- 
larged the field of quantum Hall ferromagnetism. It has 
been shora that easy-axis ferromagnetic ground states 
can occuill3 at higher filling factors when LLs with dif- 
ferent orbit r.adAw3 quantum numbers are brought close 
to alignment £3'E2l In HF language, easy-axis anisotropy 
means that many-body states with either of the two 
aligned LLs completely filled and the other empty is en- 
ergetically more favorable than the coherent superposi- 
tion state. The easy-axis pseudospin anisotropy occurs 
in this case because intra LL exchange is stronger than 
exchange between particles from LLs with different (ptbit. 
radius quantum numbers. Transport measurementsMu 
have demonstrated that easy-axis QHFs exhibit hystere- 
sis with a complicated phenomenology, presumably as- 
sociated with an interplay between disorder and domain- 
morphology similar to that in conventional thin film mag- 
nets. 

In this paper we classify QHFs according to their pseu- 
dospin anisotropy energies as either isotropic, easy-axis, 
or easy-plane systems. We report on a HF based analysis 
which predicts how the class of broken symmetry ground 
state depends on the nature of the crossing LLs. In some 
cases, competing effects allow the anisotropy energy to be 
tuned continuously generating a zero-temperature quan- 
tum phase transition between different classes of states. 
We consider only cases where no more than two LLs are 
nearly degenerate and the number of electrons is suffi- 
cient to occupy one of them. We will always assume 
that lower energy LLs, if present, are completely full 
and higher energy LLs are completely empty; coupling 
to these remote LLs can usually be treated perturba- 
tively if necessary, although we do not do so explicitly 
here. An important example of an instance in which 
more than two LLs are close to degeneracy pertains in 
double-layer systems with weak tunneling and weak Zee- 
man coupling ;L£rli3 we do not treat this or other more 
complex cases with many degenerate LLs in this paper. 
In Section II we precisely define the pseudospin language 
we use in which one of the LLs is referred to as the 
pseudospin-up state and the other LL as the pseudospin- 
down state. Since we assume that the magnetic field 
is perpendicular to the 2D electron layer, growth direc- 
tion and in-plane degrees of freedom decouple. The pseu- 
dospin quantum number then subsumes real-spin, orbit 
radius, and growth direction (subband) degrees of free- 
dom. To make the discussion more transparent we con- 
centrate on a system consisting of two nearby infinitely 
narrow 2D layers, the simplest model which has a non- 
trivial growth direction degree of freedom. Comments are 
made throughout the text about realistic samples with 
more complicated geometries. In Section III we derive a 
general expression for the HF ground state energy in the 



pseudospin ferromagnetic state. Section IV summarizes 
the rather cumbersome evaluation of Coulomb interac- 
tion matrix elements. Readers not interested in techni- 
cal details of the calculation are encouraged to skip to 
Section V where we present our conclusions concerning 
pseudospin magnetic anisotropy of single layer and bi- 
layer QHFs. This section includes phase diagrams which 
show the regimes of physically tunable parameters with 
easy-axis and easy-plane anisotropies. Symmetry break- 
ing fields and the dynamics of pseudospin reversal are dis- 
cussed in Section VI. Finally, we conclude in Section VII 
with a brief summary of the main results of our paper. 



II. PSEUDOSPIN REPRESENTATION 

The pseudospin language was introduced! to the de- 
scription of broken symmetry states in the quantum Hall 
regime, in order to draw on the analogy between double 
quantum well systems at v — 1 and 2D ferromagnets. In 
this work the pseudospin degree of freedom represented 
the layer index of a bilayer system. Here we allow the 
pseudospin index to have a more general meaning. To 
establish terminology, it is useful to recall the single- 
particle spectrum of a bilayer 2D system subject to a 
perpendicular magnetic field. Quantum well subbands 
of individual layers can be mixed by interlayer tunneling 
and shifted by the application of a bias potential. We 
limit our attention to the usual case where only the low- 
est electric subband of either quantum well is occupied 
and for explicit calculations use a zero-width quantum 
well model. Allowing for external bias and for tunnel- 
ing between the wells (see Fig. [I]), the bilayer subband 
wavefunctions of the zero-width model are 



A±i(z) 
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(1 T r A ) 1 / 2 S(z) ± (1 ± r A )^ 2 S(z - d) 



(1) 



where r A = Ay /(Ay + A 2 ) 1 / 2 , Ay is the bias poten- 
tial, A t the tunneling gap at zero bias, and d is the layer 
separation. To account for specific experimental sam- 
ples, finite width effects can be incorporated by replac- 
ing the wavefunctions (|l|) with electric subband wave- 
functions calculated using the self-consistent local-spin- 
density-approximation (LSDA) model. Ej 

In the Landau gauge the wavefunctions in the 2D plane 
take a form (j) n , s ,k(x) exp(iky)/ \JL V , where 



<t>n,s,k(x) = M 2 2 2 >!) 2 ] 1/4 H n 
(x-l 2 hf 



x - l 2 k 



x exp 



2£ 2 



(2) 



k is the wavevector label which distinguishes states 
within a LL, n — 0,1, ... is the orbit radius quantum num- 
ber, t is the magnetic length, and we have also explicitly 
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included the real-spin, s = ±1/2, degree of freedom. The 
single-particle energy spectrum consists of discrete LLs 



E. 



-|(AV 



A 2 ) 1 / 2 + htu c (r 



5) - s\g\n B B , 



(3) 

where lo c is the cyclotron frequency and the last term is 
the real-spin Zeeman coupling. Each LL has a macro- 
scopic degeneracy with the number of orbital states per 
level Ns = AB/&0, where A is the system area, B is the 
field strength, and $0 is the magnetic flux quantum. 

The many-body broken symmetry states we study in 
the following sections occur when two LLs are brought 
close to alignment while remaining sufficiently separated 
from other LLs. In our calculations, each LL can have 
one of two possible subband indices, one of two possible 
spin indices, and any value for the 2D cyclotron orbit 
kinetic energy index. The two crossing LLs can differ in 
any or all of these labels. We label one of the two levels 
as the pseudospin-up (a =|) state and the other level 
as the pseudospin down (a =1) state. We truncate the 
single-particle Hilbert space by ignoring higher LLs and 
introducing effective one-body fields that account for the 
effect of electrons in lower LLs on the two pseudospin 
states. Within this model the set of single-particle states 
reduces to following wavefunctions 

exp(iky) 



(4) 



A particle with the pseudospin oriented along a general 
unit vector rh — (sin 9 cos ip, sin 9 sin ip, cos 9) is described 
by 

9 9 

?/W-(r) = cos (-)^T,fe(^ + sin il^e^^ltkir) ■ (5) 



III. MANY-BODY HAMILTONIAN AND HF 
TOTAL ENERGY 

In the HF approximation, the QHF has a single Slater 
determinant state with the same pseudospin orientation 
for every orbital k. In this section we derive general ex- 
pressions for the dependence of the many-electron state 
energy on pseudospin orientation. It is convenient to 
express the many-body Hamiltonian using Pauli spin- 
matrices t x , T y , and t z and the 2x2 identity matrix 
which we label T\. In this representation the Hamilto- 
nian reads 



H 



iv 2 

E E E hiT i ' ot( i(c t '),k c <o l ),k 



i=X,x,y,z k—1 a,a f = l 



E 



k 2 ,k> =1 



2 

E,„k'- L ,k' 2 ,k\,k 2 a^Qi a' 2 -a2 



1' 



X C i(a[),k[ C l(a' 2 ),k' 2 C c'{a2)M C <?(ai),k 1 , 



(6) 



where cr(l) =T and cr(2) =J.. The one-body terms 6, 
include, in general, the external bias potential, tunneling, 
cyclotron and Zeeman energies and also the mean-fields 
from interactions with electrons in the frozen LLs lying 
below the a —\ and J, levels. We will give an explicit 
expression for bi in section VI. Here and in the following 
two sections, we concentrate on the two-body terms in 
the Hamiltonian @. 

The potentials Wij represent different combinations of 
Coulomb interaction matrix elements, V„> „> ai ^ 2 , of the 
single-particle pseudospin states 

e 2 

x -p ^V'ffi.fci^OVVa.kaOrb) (7) 

General expressions for the pseudospin dependent inter- 
actions Wij are given in Table | in terms of the following 
matrix element combinations: 



Bf = 


\^ 


,T,T.T 


±n,Lw) 


,Bf = 


T^T.LTd 


± VI,t,j.,t) 


Bt = 


1^ 


,W>T 


± v i,r,r,i) 


,B± = 


1 


± VLl,t,t) 


Bf = 


\^ 


,T,T,1 


± v i,i,i,i) 


,Bt = 


1 


± v l,l,l,1) 


B± = 


\M 


,T,Lt 


± VI,t.t,t) 


,Bt = 


1 

7(^,1, 1,1 





(8) 

In (||) we have omitted orbital guiding center indices for 
simplicity. 

The many-electron state with pseudospin orienta- 
tion to is |*[m]) = l\k=i4h,k\°}' where c A,fc creates 
the single-particle state whose wavefunction is given in 
Eq. (§). We find that 



eHFifh) = 



{^[m]\H\^[rh]) 

n2 



w, 



i—x,y,z 

2 E ■ '■' 



where 



1 



9 fel,fe 2 =l 



( U/ ,k i, k 2,k 1 ,k 2 



k 2 ,kx,ki,k 2 



(9) 



(10) 



has direct and exchange contributions. Eq. (Q) is the 
most general form for the HF energy of a pseudospin-1/2 
QHF. The magnetic anisotropy of the particular QHF 
system is governed by the terms in @ that are quadratic 
in the pseudospin magnetization rrii. The values of the 
coefficients U^j depend on the nature of the crossing LLs 
and and can result in isotropic, easy-plane, or easy-axis 
quantum Hall ferromagnctism. 
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IV. PSEUDOSPIN MATRIX ELEMENTS OF THE 
COULOMB INTERACTION 

This section contains the derivation of explicit expres- 
sions for the anisotropy energy coefficients Uij, assuming 
the zero- width quantum well model wavefunctions (0). 
Using the Fourier representation of the Coulomb interac- 
tion we write the pseudospin matrix elements as 



,rk' 1 ,k' 2 ,ki,k 2 _ c 



a / f "iiy,K 1 —K 1 "—q v ,k 2 —k2 

q 

x e lq ^ k i +kl)/2 e- iq * (k i +k2 *> /2 v a[ y 2 , aua2 {q) (11) 
and hence 

- \ t-trki,k 2 ,ki,k 2 yk 2 ,ki,ki,l 

/ j \ <j' 1 ,<7 2 ,a 1 ,a 2 tr[, cr' 2 , rrx,o 



,0"2 / 



Aii ,k 2 — 1 



v^a'^u^iO) - -jY^v^^^iq) (12) 



In the following we take £ as the unit of length and e 2 /el 
as the unit of energy. Then N$/A = l/2it and Eq. jl^ ) 
can be rewritten in the more transparent form 



_ \ °~ (yki,k 2 ,k u k 2 _ yk 2 ,ki,ki,k 2 \ 



9 fci,fe 2 =i 

T^j^e-^^/^^^CQ)-^,^^,^^] . (13) 

The first factor in square brackets in this equation origi- 
nates from the Hartree contribution to the energy while 
the second factor originates from the exchange contribu- 
tion. 

In these equations v a ^ a ^ ! j 1 ^ 2 (q) is a pseudospin- 
dependent effective 2D interaction which, because of the 
separability of the in-plane and out-of-plane degree-of- 
freedom terms in the single electron Schroedinger equa- 
tion, is the product of two factors: the subband factor 



J — CO J — CO 

x A «K)( z i) A «K)<>2) A ?((Tl) (zi) A ?((T2) (z 2 ) (14) 
and the in-plane term 



■j n , , 

CT, ,fT 2 ,CTl ,<?2 v 



,(*) = e? 

X / da:i0n(tri).*«).?v/2( :El )^n(ffi).8(ffi).-g»/2( a; l) 
</ — oo 

/"CO 

X / d,X 2 <f)n(a' 2 ),s(<?'2) -qy/2(x2)K(a 2 ),s(a 2 ),q y /2{x2) (15) 



For a general sample geometry the subband factor 
( |j"4| ) has to be calculated numerically using the self- 
consistent LSDA wavefunctions. For our model bilayer 



system, however, we can obtain analytic expressions for 

V a> a'^au^iO)- 111 the CaSe ° f ^( CT ) = £( _cr )> E q S - (S) and 

@ give 



v%^a 2 =\[^+rl) + {l-rl)e- dq ] . (16) 



In the second case, i.e. £(er) = — £(— cr) 

1 

2 



^, CT . CT , CT -7T[(l + 4) + (l-r|) e - d9 ] 



<- CT)< r,- CT = 5 K 1 - r i) + (i + Ay 



dq] 



and 



v%^ 2 =v r -j(l-rl)^(l-e- dq ) 

if + N =±i. 



(17) 



For the in-plane factor in the effective interaction it is 
necessary to distinguish several cases. If the pseudospin- 
up and pseudospin-down levels have the same real- 
spin index and same orbit-radius quantum number, i.e. 
s(<r) = s(— cr) and n(a) = n{—cj) = n, the effective inter- 
action is independent of the pseudospin indices and we 
obtain from (g) and ( |l5| ) 



N 

V rr' 1 ,a' 2 ,cr 1 .cr 2 



27T 



(18) 



where L n (x) is the Laguerre polynomial. For identi- 
cal spins but different orbit-radius quantum numbers, 
i.e. s(a) = s(—a) and n(a) ^ n{—a), we define n < = 
min [n(cr), n(— cr)] and n> = max [n(cr), rt(-cr)] . Then 
this factor in the effective interaction is 



.•IV 



.N 



2tt 



L n(a) (q 2 /2)] 2 , 



2tt 



fj, — cr, — (T,CT 



27T ro<! / q 



I /„2\™>-"< 



9 n>! 



! V 2 



[^>-<(g 2 /2)] ; 



otherwise « CT ^ >CTl)CT2 = 



(19) 



If the two pseudospin LLs have opposite real-spins then 
only scattering processes that conserve pseudospin (and 
therefore also real-spin) at each vertex will contribute to 
the anisotropy energy. For s(cr) = — s(— cr) and n(cr) = 
rt(— cr) = n we obtain 



a. a, a, a ^a. — a.a. — a 



otherwise v^y^^^ 



2?r 

q 

= o 



in(g 2 /2)] 2 , 



(20) 
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Finally if the pseudospin LLs have opposite real-spins 
and different orbit-radius quantum numbers, the in-plane 
factor in the effective interactions is 



a, a. a, a 



.N 



y[W<r72)r 

2tt 



a, — cr.cj, — a 



otherwise 



— L n{a) {q 2 /2)L n{ _ a) {q 2 l2) 



T 2 ,{Ti,CT2 



(21) 



for the subband factor and explicit expres- 
(|1S|)-(21) for the in-plane factor in the effective in- 
teraction, together with Eqs. @,@,(|) and Table @, 
provide a formal recipe to calculate the anisotropy co- 
efficients Ui.j for crossing LLs with any combination of 
quantum well subband, orbit radius and real-spin inde- 
ces. In the following section we use the explicit forms 
( 17]) and ( |l6| ) for the subband factor to develope a pseu- 
dospin anisotropy classification scheme for our model bi- 
layer system. 



V. MAGNETIC ANISOTROPY 

The nature of the anisotropy energy is of qualitative 
importance for two-dimensional ferromagnets, including 
QHFs. Systems with easy-axis anisotropy, i.e., discrete 
directions at which the energy of the ordered state is min- 
imized, have long range order at finite temperature and 
phase transitions in the Ising universality class. Systems 
with easy-plane anisotropy, i.e., a continuum of copla- 
nar pseudospin magnetization orientations at which the 
energy of the ordered state is minimized, do not have 
long-range order but do have Kosterlitz-Thouless phase 
transitions at a finite temperature. In the isotropic case, 
all directions of pseudospin magnetization have identical 
energy, only the ground state has long-range order, and 
there are no finite-temperature phase transitions. Nev- 
ertheless, magnetization correlations become extremely 
long at low temperatures. The objective of this work is 
to predict which class of QHF occurs for a particular pair 
of crossing LLs. 

We start our analysis by considering pseudospin LLs 
that belong to the same subband, i.e. £(t) = £(!)■ In this 
case, only LLs with opposite real-spins can be aligned. 
Two examples of QHFs falling into this category, total 
filling factor v = 1 with n(|) = n(l) = 0, and v = 2 with 
n(t) = 1, n(l) = 0, are illustrated in Fig. ||. The cases 
of 7i(T) 7^ "(I) are realized when the ratio of the spin- 
splitting to LL separation is an integer. In GaAs, this 
ratio is only ~ 1/60 at perpendicular fields but can be 
tuned by tilting the magnetic field away from the normal 
to the 2D layer. For typical well widths orbital effects of 
the in-plane field, not included here, become important 
at the tilt angles where the coincidences of interest are 
realized and have to be accounted forlifj to obtain correct 
values of the pseudospin anisotropy energy coefficients 



(|Io|). However, recent work on AlAs quantum wellsS 
and InSb quantum wellsEa with large Zeeman couplings 
have made the situation we study below, which assumes 
perpendicular magnetic field, accessible. 

When opposite spin LLs cross, Eqs. ( pfj| ) and (|2l]) im- 
ply that all interactions Bf with i > 2 in (||) vanish. 
Then the only non-zero anisotropy term is 

1 f 00 

0*,* = -8 J d q e-i 2 / 2 [L nm ( q 2 /2)-L n(L) ( q 2 /2)] 2 
x [(l + rl) + (l-ri)e-*] . (22) 

Note that the Hartree energy contribution to anisotropy 
always vanishes when the crossing LLs share the same 
subband wavefunction. If the two pseudospin levels 
also have the same orbit-radius quantum number then 
Eq. ( p2| ) gives U z ,z = and the ferromagnetic state is 
isotropic. Physically, the result follows from the indepen- 
dence of the Coulomb interaction strength on real-spin. 
An important example of these isotropic QHFs occurs 
when n(T) = = and rj\ = 1, i.e., there is no tunnel- 
ing between layers-This is the thoroughly studied single- 
layer v = 1 QHFI3 □ for which the HF theory ground 
state happens to be exact. We remark that quantitative 
estimates based on the HF mean-field theory presented 
here require corrections to quantum fluctuation effects in 
cases when the ordered pseudospin moment direction is 
not a good quantum number. A detail understanding of 
these corrections is one challenge for future experimental 
and theoretical work on QHFs. 

For n{\) ^ n(i), Eq. (§|) implies that U z , z < 0, 
making the z-axis the easy pseudospin orientation axis. 
Again, at Ta = 1 our model reduces to that of a sin- 
gle layer 2D systems whose easy- axis anisotcopy at even 
filling factors has been identified previously]!!] In finite- 
thickness single quantum wells, a QHF with pseudospin 
LLs of the same subband but different real-spin and orbit- 
radius indices is also easy-axis. The magnitude of the 
anisotropy will decrease with layer thickness, as can be 
seen by comparing U z>z in (|2|) calculated for ta. = 1 and 
ta = 0. (Note that the single-subband unbiased double 
well with finite tunneling, i.e. ta = 0, models a single 
layer system with an effective thickness d.) 

We now turn to the crossing of LLs with different 
subband indices, for which the pseudospin anisotropy 
physics is richer. In Fig. |^ we show examples of n(\ 
) = n (l) bilayer QHFs for v = 1 and v = 2 based on 
same real-spin and opposite real-spin LLs respectively. 
Eqs. @,@,(|j^,(jl!), and Table | imply four non-zero 
anisotropy terms for n(|) = n(|) = n and s(|) = s(\): 



'A i 



V x , x = u(l - r 2 A ) , 



where 
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U= 2-~2 



dge-9 a / a [L n (9 2 /2)] 2 (l-e- d «) . (23) 



First term in the expression for energy u comes from the 
Hartree interaction, the second term represents exchange 
contribution which is always smaller than the Hartree 
energy in this case, i.e., u > 0. For zero tunneling be- 
tween layers (r A = 1), the only non-zero anisotropy en- 
ergy component, U zz — u, is positive leading to the easy- 
plane anisotropy of the QHF. In the absence of symme- 
try breaking fields (the linear pseudospin magnetization 
terms in (|9|)) the variational energy (|^) is minimized when 
the pseudospins condense into a state magnetized at an 
arbitrary orientation within in the x — y plane. In this 
state, electronic charge is distributed equally between the 
layers (pseudospin angle 9 = 0) minimizing the electro^ 
static energy; spontaneous interlayer phase coherences 
the physical counter part of pseudospin order in this case, 
lowers the total energy of the system by strengthening in- 
terlayer exchange interactions. 



For r A < 1, Eqs. (£23|) imply the following quadratic 
terms in the HF total energy 

Uijiriinij = u[r A m z + (1 - r 2 A ) 1/2 m x ] 2 , (24) 



i.e., the easy-plane is tilted from the x — y plane in the 
pseudospin space by angle a = arctan [(1 — r^) 1 / 2 /Va] • 
The pseudospin basis states at different values of r A are 
related, however, by a unitary transformation, which cor- 
responds precisely to a rotation about y-axis by angle —a, 
as seen from Eq. ([!]). The easy-plane where the above 
anisotropy energy is constant is, for any value of r&, the 
plane of equal charge per layer. 

When pseudospin-up and pseudospin-down states dif- 
fer by more than their subband indices, by their spin 
indices for example, the two pseudospin basis sates at 
different r A are not related by a unitary transforma- 
tion. In this case the magnetic anisotropy does depend 
on r A . For example, consider the case f(f) = — 
n(t) = n(i), and s(T) = — sQ). For opposite real-spin 
LLs, all anisotropy terms that include pseudospin non- 
conserving scattering processes drop out. The only non- 
zero energy term, U zz , has the same value as in (23). 
Hence, the QHF is isotropic in the unbiased (ta = 0) bi- 
layer system while applying external bias (r& > 0) leads 
to easy-plane anisotropy in pseudospin space. 

At this point let us make an experimentally important 
comment on thej-,biLayer systems realized in wide single 
quantum wells .Hut^l The difference between this sample 
geometry and the double quantum well with narrow (in 
our model infinitely narrow) layers is in the nature of 
the barrier responsible for the bilayer character of the 
electronic system. In wide quantum wells the barrier 
is softjEj originating from Coulomb interactions among 
electrons in the well. Then the tunneling probability be- 
tween layers is strongly dependent on the electron density 
and quantum well subband populations. This tunabil- 
ity makes wide single quantum wells an experimentally 



attractive alternative to double wells in studies of bi- 
layer quantum Hall phenomena. The softness of the bar- 
rier can, however, lead to qualitatively important conse- 
quences for the ordered many-particle states. Translated 
into the pseudospin language, A t cannot be treated as an 
external one-body field acting-pn the pseudospin particles 
but, in general, will depencO on the pseudospin orien- 
tation in the ordered ground state. For the pseudospin 
LLs discussed in the previous paragraph (£(T) = — 
n(|) = n(l), and s(f) = — s(|)) and for r A — 0, this ef- 
fect can lead to an anisotropic QHF. LSDA calculations 
indicate that at low electron densities the anisotropy will 
be easy-plane while easy-jaxis anisotropy is more likely to 
develop at high densities O We make this remark to point 
out that not all results obtained for the double quantum 
well model are directly applicable to bilayers in wide sin- 
gle wells. In many cases the theoretical description of 
QHFs in wide single wells requires modifications of the 
idealized bilayer model to account for mixing of higher 
electrical subbands. The self-consistent LSDA for the 
growth direction single-particle orbitals is a particularly 
convenient, if somewhat ad hoc, method that allows any 
sample geometry to be studied while retaining the basic 
structure of the many-body HF formalism for QHFs. 

In the remaining part of this section we consider pseu- 
dospin LLs with opposite subband indices and differ- 
ent orbit radius quantum numbers. At total filling fac- 
tor v = 3, for example, the pseudospin LLs will have 
the same real-spin while at v = 4 opposite spin LLs 
can be aligned, as shown in Fig. |ij A common feature 
of the QHFs discussed below is the transition from a 
state with easy-axis anisotropy to a state with easy-plane 
anisotropy as r A and the layer separation d are varied. 
For s(t) = s(l), Eqs. © , @ , @ , © , and Table § give 
three non-zero anisotropy energies 



^ -I «v 



-q 2 /2 



x|[L„ (T) (g 2 /2)-i n(i) ( g 2 /2)] 2 (l + e- d «) 
+r 2 A [L n{]) (q 2 /2) + L n(i) {q 2 /2)}\l - e^)} , 

2 \ n^-n. 



g2/2 n<! / <r 
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U x>x = U v ,y = -- I dqe 

x[L->- n <{q 2 /2r[\ - dsn ' ••") . 

where 



n< = mm 



[n(T),n(l)]> n>=max[n(T),nU)] . (25) 

The HF total energy contributions that are quadratic 
in the pseudospin magnetization components can be 
grouped as 



Ui^rmmj = (U Z}Z ~ U XtX )ml + U XjX . (26) 



G 



(Recall that rh is a unit vector, i.e. m 2 + m 2 = 1 — m 2 z .) 
From Eq. @ we obtain that for U z . z -U x . x < the QHF 
has easy-axis anisotropy while for U ZiZ — U x<x > the 
system is an easy-plane ferromagnet. At the critical layer 
separation d = d* , obtained from the condition U Z}Z — 
U XjX , the magnetic anisotropy vanishes and a fine-tuned 
isotropy is achieved. It follows from Eqs. ( p5| ) that d* is 
finite for all values of r A . 

For pseudospin LLs with s(|) = — the anisotropy 
energy components U XiX and U y ^ v vanish and the critical 
layer separation d* corresponds to U xx = 0, where U z<z 
is given by the same expression as in the s(|) = s(i) case 
(see Eq. (p5|)). Since U ZjZ < at r A — 0, the critical 
separation d* diverges in the absence of external bias, 
i.e., easy-axis pseudospin anisotropy does not exist for 
any layer separation. For r A > the transition between 
easy-plane and easy-axis anisotropy occurs at finite d as 
for the s(t) = s{{) pseudospin LLs. 

In Figs. g|(a) and|5|(b) we show the magnetic anisotropy 
phase diagrams in the d-r\ plane calculated for v = 3 and 
v = 4 QHFs (see Fig. Since the layer separation is in 
units of magnetic length, these figures imply that transi- 
tion between easy-axis and easy-plane anisotropics at a 
given filling factor can be induced in one physical sample 
by changing the density of the 2D electron system. High 
electron densities would correspond to the easy-plane re- 
gion, and low densities to the easy-axis region. Note that 
these numerical results confirm the general remark made 
above, since the critical layer diverges as r A — > for 
v = 4 while it remains finite for v = 3. 



VI. SYMMETRY BREAKING FIELDS 

The pseudospin orientation in a QHF ground state 
is determined by minimizing the variational total en- 
ergy (^). In the absence of energy terms that are lin- 
ear in pseudospin magnetization components, the HF or- 
dered states spontaneously break continuous SU(2) or 
U(l) symmetry^ in the case of isotropic or easy-plane 
QHFs respectively, and the discrete symmetry between 
pseudospin-up and pseudospin-down orientations in the 
case of easy-axis QHFs. In this section we take into ac- 
count external and internal potentials which contribute 
to the linear terms in the HF total energy and comment 
on the pseudospin reversal that can be triggered by ad- 
justing these symmetry breaking fields. We focus on an 
case which we feel is particularly appropriate for experi- 
mental study by considering the ordered v = 3 quantum 
Hall state. Similar considerations would apply for all 
classes of QHFs discussed in this paper. 

The pseudospin LLs in the bilayer v = 3 QHF (see 
Fig. ^) have opposite subband indices and orbit radius 
quantum numbers n = and n = 1, respectively. We 
call the [£ = — l,n = 0,s = +1/2] LL the pseudospin-up 
state and the [£ = l,n = l,s = +1/2] LL the pseudospin- 
down state. With this definition, the one-body potentials 



in (|9|) can be written as 

b z = -~[(A 2 r + A 2 ) 1 / 2 -?uJ c + I F -I H , z ] , 

b v =Q. (27) 

The effective field If is the difference between 
pseudospin-up and pseudospin-down particle exchange 
energy with electrons in the fully occupied [£ = l,n = 
0,s = +1/2] LL, i.e, 




For ?'a > 0, electrons in the [£ = l,n = 0,s = ±1/2] 
LLs produce also an electrostatic field, represented by 
Ih,z and Ih,x in (p^), which screens the external bias 
potential. Since this effective field favors occupation of 
a particular layer rather than a particular pseudospin 
state it couples to z and x components of the pseudospin 
operator. The energy inbalance between the two layers 
produced by Ih is 2drA which, together with Eq. (Q), 
gives 

I H , Z = 2dr 2 A , 

Ih,x = 2dr A (l-r 2 A )^ 2 . (29) 

In the expression for the HF total energy (^), we in- 
cluded explicitly the contribution to the symmetry break- 
ing fields which results from Coulomb interactions be- 
tween electrons in the pseudospin LLs. For the v = 3 
QHF we are considering here, only U% jZ and U Zt \ ener- 
gies are non-zero: 

U ltZ = U z ,i - -- / dqe-fl 2 [1 + r\ + (1 - rl)e- d «] 

x (q 2 - 9 4 /4) • (30) 

In the bilayer system with no tunneling (r A = 1), 
Ih,x — and the total symmetry breaking field, b* = 
b z — Ui tZ , is oriented along the z pseudospin direction 
and is given by 

b* = -- [(A^ + A 2 ) 1 / 2 - Tiu c + V^72/2 - 2d 
- ■ (31) 

In Figs. ||(a)-(c) we plot the pseudospin evolution with 
effective field b* for the three anisotropy regimes of a 
v = 3 QHF with r A = 1. At the phase boundary be- 
tween easy-axis and easy-plane anisotropies, the v = 3 
QHF is isotropic and the pseudospin reverses abruptly 
at b* — (see Fig. ||(a)). In the easy-plane anisotropy 
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regime, (U z<z — U XiX > 0), the pseudospin evolves contin- 
uously with b* as illustrated in Fig. §(b) reaching align- 
ment with b* at |6*| > U z<z — U x . x . For the easy-axis 
anisotropy case (U z z — U x _ x < 0), the HF energy has two 
local minima at m z = ±1 when < \U Z>Z — U XjX \. The 
pseudospin-up and pscudospin-down polarized states are 
separated by an energy barrier which results in the hys- 
teretic pseudospin-reversal behavior shown in Fig. [|(c). 

In bilayer systems with non-zero tunneling, pseudospin 
reversal follows a more complicated pattern in which the 
competition between x and z components of the sym- 
metry breaking field plays an important role. In gen- 
eral, the pseudospin will rotate in the x-z plane, i.e. 
m x = yl — m 2 z . Since the derivative of the HF energy 
with respect to m z diverges at m z = ±1 due to the Ih,x 
term the pseudospin will never align completely with the 
z-axis when ta < 1. 

VII. SUMMARY 

In the strong magnetic field limit, the physics of 
high mobility two-dimensional electron systems is usu- 
ally dominated by electron-electron interactions except 
at integer filling factors, where the single-particle physics 
responsible for the gap between Landau levels assumes 
the dominant role. When external parameters are ad- 
justed so that two or more Landau levels simultaneously 
approach the chemical potential, the integer filling fac- 
tor case is less exceptional, interaction effects are always 
strong, and uniform density broken symmetry ground 
states analogous to those in conventional ferromagnets 
are common. In this paper we have discussed how the 
nature of these states depends on the character of the 
nearly degenerate Landau levels. Our attention is re- 
stricted to the case where only two Landau levels are 
close to the chemical potential and we distinquish these 
crossing Landau levels by introducing a pseudospin de- 
gree of freedom. 

Using Hartree-Fock variational wavefunctions, we are 
able to derive an explicit expression for the dependence of 
ground state energy on pseudospin orientation for cross- 
ing LLs with any combination of quantum well subband, 
orbit radius and real-spin degree-of-frcedom quantum 
numbers. As in conventional magnetic systems, quali- 
tative differences exist between the physical properties 
of isotropic (Heisenberg) systems with no dependence 
of energy on pseudospin orientation, easy-axis (Ising) 
systems with discrete prefereed pseudospin orientations, 
and (XY) easy-plane systems for which the minimum is 
achieved simultaneously for a plane of orientations. Our 
mean- field results predict which class of pseudospin quan- 
tum Hall ferromagnet occurs in different circumstances. 
We focus on a model commonly used for bilayer quantum 
Hall systems in which the finite width of both quantum 
wells is neglected. The external parameters of the model 
arc the Zccman coupling strength, the bias potential be- 



tween the wells Ay, and the single particle splitting due 
to interlayer tunneling A t . In the limit A t = 0, this 
model applies to a single quantum well when the cross- 
ing Landau levels have the same subband wavefunction, 
i.e., are in the same quantum well. Classification predic- 
tions for this model as a function of layer separation d 
and the ratio ta = Ay/(A^ + Af) 1 / 2 are summarized in 
the diagram |^. 

For the single quantum well case we find isotropic be- 
havior when the orbit radius quantum numbers of the 
crossing Landau levels are identical, and easy-axis behav- 
ior otherwise. In general when the crossing Landau levels 
have identical subband wavefunctions, the nature of the 
pseudospin anisotropy does not depend on the paramters 
d and ta • For different subband wavefunctions the pseu- 
dospin anisotropy can vary in the d— ta plane. Particu- 
larly intriqueing is the case of crossing LLs with different 
subband and orbit radius quantum numbers where, at a 
given filling factor, the system can undergo a quantum 
phase transition from an easy-axis to easy-plane QHF. At 
the phase boundary the pseudospin anisotropy vanishes 
and a fine-tuned isotropy is achieved. The critical values 
of d and ta are experimentally accesible and may be ac- 
companied by observable changes in pseudospin reversal 
properties as external parameters are varied. 
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TABLE I. Coulomb interaction matrix elements Wi, 
i,j — 1, x, y, z. The B^ terms are defined in (H). 



FIG. 1. Schematic of the conduction band edge profile of 
a biased double-quantum-well sample with non-zero tunnel- 
ing between 2D layers. Energy levels due to quantization of 
electron motion along the growth direction are also indicated. 
For concrete calculations we assume infinitely narrow quan- 
tum wells separated by the distance d. 



FIG. 2. Schematic LL diagrams for v = 1 (a) and v — 2 
(b) single-subband QHFs. Index o =f , J, labels the pseudospin 
LLs. The crossing LLs are indicated by half-solid dots, while 
inert filled LLs are indicated by solid dots and inert empty 
LLs by open dots. 



FIG. 3. Schematic LL diagrams for v — 1 (a) and v = 2 
(b) bilayer QHFs. The pseudospin LLs have opposite subband 
indices and same (a) or opposite (b) real-spins. 



FIG. 4. Schematic LL diagrams for v — 3 (a) and v = 4 
(b) bilayer QHFs. The pseudospin LLs have opposite subband 
indices, different orbit radius quantum numbers and same (a) 
or opposite (b) real-spins. 



FIG. 5. Magnetic anisotropy phase diagrams for bilayer 
v = 3 (a) and v = 4 (b) QHFs from Fig. §. In the white 
region the anisotropy is easy-axis, in the grey region the QHF 
has easy-plane anisotropy. At the phase boundary the ferro- 
magnetic state is isotropic. 



FIG. 6. Pseudospin orientation as a function of the effective 
field b* for the isotropic QHF and as a function of b* relative 
to the anisotropy energy \U Z)Z — U x ,x\ for the easy-plane (b), 
and easy-axis (c) QHFs at v = 3. 



FIG. 7. Magnetic anisotropy of QHFs with pseudospin LL 
subband indices £(cr), orbit radius quantum numbers n(er), 
and real-spins s(o). 
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